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> ' ABSTRACT 

o 



We study r-modes of a rotating magnetized neutron star, assuming a magnetic dipole 
field whose axis is aligned with the axis of rotation. We approach the problem by 
' applying a singular perturbation theory to the oscillations of rotating stars. In this 

CO \ treatment, we divide the star into a thin surface magnetic layer and a non-magnetic 

core. We integrate linearized ideal MHD equations in the surface magnetic layer and 
non-magnetic oscillation equations in the core, and match the two integrations at the 
interface to obtain a complete solution. For a polytropic neutron star model of mass 
\ M = IAMq and radius R — lO^cm, the magnetic dipole field becomes effective on 

' the modal properties of the r-modes only when the field strength Bs is much greater 

than lO^^G. We also find that the damping effects caused by very short magnetic 
■ ■ perturbations in the surface layer are not important for the r-mode instability of 

\ rapidly rotating neutron stars if the field strength Bs is smaller than lO^^G. 

f~i ■ Key words: instabilities - stars: neutron - stars: oscillations - stars : rotation - stars 

C^l : magnetic fields 
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1 INTRODUCTION 
> 

• ^ • Several numerical investigations have recently been carried out to study the properties of acoustic modes of magnetized 
' stars (e.g., Dziembowski & Goode 1996, Bigot et al 2000, Cunha & Gough 2000, Saio & Gautschy 2004), particularly for 
^ . roAp stars (Kurtz 1990). These authors solved the problem of the oscillations of a magnetized star by applying a singular 
■ ■ ' perturbation theory to the stellar oscillation. In their treatment, the stellar interior was divided into the thin surface magnetic 
layer where the magnetic pressure pm. = -B^/Svr is comparable to or greater than the gas pressure p and the non- magnetic 
core where the gas pressure is dominant over the magnetic pressure (Biront et al 1982, Roberts & Soward 1983, Cambell 
& Papaloizou 1986). Linearized ideal MHD equations were integrated in the thin surface magnetic layer, and non-magnetic 
oscillation equations were integrated in the core, and then the two integrations were matched at the interface between the 
magnetic and non-magnetic regions to make a complete solution. This treatment is based on the following properties of 
magnetic perturbations in the interior of stars. If we assume that the gass pressure dominates the magnetic pressure, or 



equivalently that the Alfven wave velocity va = \/ B^/inp is much slower than the sound velocity vs ~ \/p/p in the deep 
interior, it is reasonable to consider that the magnetic perturbations are largely determined by the induction equation. But, 
there also exists a component of magnetic perturbations that originates from the Lorentz force term j x B in the equation of 
motion. For a typical oscillation frequency of acoustic modes, this component appears as very short magnetic perturbations, 
the wavelengths of which are much shorter than those of the acoustic modes. In the limit of /3 = p/pm ^ oo in the deep 
interior, we can expect the short wave component of the magnetic perturbations may be lost from the wave system suffering 
strong dissipation, which makes it legitimate to divide the stellar interior into a magnetic surface layer and a non-magnetic 
core and to disregard the short wave component in the core (e.g., Biront et al 1982, Roberts & Soward 1983). If we assume 
that the oscillation energy imparted by the short magnetic perturbations is lost from the wave system (Biront et al 1982, 
Roberts & Soward 1983), this loss of the magnetic perturbation energies works as a damping mechanism for the acoustic 
modes. This damping mechanism was implemented in global acoustic mode calculations, for example, by Bigot et al (2000), 
Cunha & Gough (2000), Saio & Gautschy (2004). 
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Although the studies mentioned above were all for non-rotating stars, Shibahashi & Takata (1993) and Takata & Shiba- 
hashi (1995) have developed a theory of the oscillations of rotating magnetized stars, treating both rotation and the magnetic 
field as small perturbations to the oscillations of non-rotating and non-magnetic stars. They recognized, however, that the 
perturbational treatment of the magnetic field in the surface layer might not be justified. Bigot & Dziembowski (2002) have 
also considered the effects of a magnetic field and rotation on acoustic oscillations. Although they assumed a slow rotation to 
treat high radial order p-modcs, they employed a non-perturbative approach for the magnetic field, to improve the analyses 
by Shibahashi & Takata (1993) and Dziembowski & Goode (1996). 

J?-modes of neutron stars have attracted good attention recently because of the recognition that the r-modes are driven 
unstable by the emission of gravitational radiation (Andcrsson 1998, Friedman & Morsink 1998) duo to Chandrasekhar- 
Friedman-Schutz instability mechanism (Chandrasekhar 1970, Friedman & Schutz 1978 a,b). The r-modes of rotating neutron 
stars are expected to be a good candidate for a gravitational wave source (e.g., Andersson & Kokkotas 2004). They are also 
regarded as a mechanism that decelerates the spin of young neutron stars, since the gravitational radiation carries away 
angular momentum of the star (e.g., Owen et al 1998). However, before we draw a decisive conclusion concerning the role 
played by the r-modes in neutron star physics, we need to examine how strongly the modal properties are affected by various 
physical factors such as viscosity, solid crust, superfluidity, general relativity, nonlinear amplitude saturation, and so on. A 
magnetic field of the star is one of such influential factors. 

i?-mode oscillations of rotating magnetized neutron stars have been discussed by RezzoUa, Lamb, & Shapiro (2000), Ho & 
Lai (2000), and Morsink & Rezania (2002). For example, Ho & Lai (2000) estimated the shifts of the r-modc frequency caused 
by the strong magnetic fleld using a local linear analysis, suggesting that the oscillation manifests itself as Rossby waves in the 
bulk interior, but assumes the Alfven wave character near the stellar surface. Morsink & Rezania (2002) developed, based on 
the formalism introduced by Schenk et al (2002), a formalism to compute the magnetically modified global r-modes, without 
assuming a weak magnetic field of the star. In this paper, we approach the problem of the oscillations of a rotating magnetized 
star by applying a singular perturbation theory to the oscillations of a rotating non- magnetic stax (e.g., Lee & Saio 1986). In 
§2, we derive the oscillation equations of a rotating magnetized star, assuming the axis of the dipole field is aligned with the 
axis of rotation. §3 is for numerical results, and §4 is for discussions and conclusions. 



2 METHOD OF SOLUTION 

2.1 Oscillation equations of a rotating star with a magnetic dipole field 

In this paper, we assume adiabatic oscillations of a uniformly rotating star and employ the Cowling approximation, neglecting 
the Euler perturbation of the gravitational potential. We also assume that the axis of the magnetic dipole field is aligned with 
that of rotation. If we employ spherical polar coordinates (r, ^, 0) whose origin is at the center of the star, a magnetic dipole 
field may be given by 

B = /iV (cos6>/r^) , (1) 

where ji is the magnetic dipole moment, and we have assumed that the axis of the field corresponds to the coordinate axis 
defined by 6 = 0. When the ax;es of rotation and the magnetic field are aligned, it is possible to assume that the equilibrium 
state is axisymmctric with respect to the axe = Q, and that the time and azimuthal angular dependence of the perturbations 
is given by ^ where cr is the oscillation frequency in an inertial frame and m is an integer representing the azimuthal 

wave number. Note that the magnetic dipole field is a force free field (i.e., (V x B) x B = 0) and does not affect the equilibrium 
structure of the star. Since we are interested in low frequency oscillations in this paper, we neglect the effects of rotational 
deformation. Linearized ideal MHD equations with infinite conductivity may be given by (e.g., Unno et al 1989) 



- o)^^ = — Vp' + -hrVp - 2iujn — (V X S') X -B, (2) 

p' + V • (p4) = 0, (3) 

B' = V X (^ X B), (4) 

and for adiabatic pulsations 

where p is the mass density, p is the pressure, fl is the angular velocity of rotation, B is the magnetic field, and the physical 

quantities with a prime (') denote their Euler perturbations, ^ is the displacement vector, and uj — a + mCl denotes the 
oscillation frequency observed in the corotating frame of the star. Here, A is the Schwartzshild discriminant defined by 

. _ dln/9 1 dlnp , , 
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where 



Ti = 



dlnp 
d\np I 

' / ad 



(7) 



To represent the angular dependence of the perturbations of a rotating magnetized star, we expand the perturbed 
quantities in terms of spherical harmonic functions with different degrees I for a given m. The displacement vector ^ is then 
given by 



oo r 



and the perturbed magnetic field B' is given by 

oo 

B', = Bo{r)y^bf,{r)Y{r{0,cl>)e"^\ 



OC r 



S; = So(r)^ 



and a perturbed scalar quantity such as p' is given by 

oo 

p' = J2pi3(^)Y{;{e,4>)e^''\ 



(8) 
(9) 
(10) 

(11) 

(12) 
(13) 

(14) 



where Bo{r) = n/r^ , and Ij = \m\ + 2{j — 1) and I'j = Ij + \ for even modes, and Ij = \m\ + 2j — \ and I'j = Ij — 1 for odd 
modes, respectively, and j = 1, 2, 3, ■ ■ ■. 

Substituting the expansions into the linearized basic equations (2) to (5), we obtain a set of coupled ordinary linear 
differential equations of infinite dimension for the expansion coefficients. If we use for the expansion coefficients vector notation 
given by 

{Vih=SiJr, {y,)j=pypgr, {b%=bf,,, {b")j = bP,, {b^)^=hl, {h)j = Hi./r, {t)j=Ti,/r, (15) 
the oscillation equations of a rotating star with an aligned magnetic dipole field are given in the Cowling approximation as 
= {Vg - 3) yi - Vay^ + ^oh, (16) 



r-j^ = (ciw^ + rAj j/j + (1 — rA — U)y2 — 2mciu)Qh — 2ciujflCoit 



dr 



^ db" ^ dib-r _ b" Ab-r ^ b^ 

Cor- h mr- /Go mf Go — 

dr a dr a a a a 



(17) 



d ( b" /a \ _ ( miyG:[\K + Mo) - AiK/a G-^MoAo/ci 



\ I Vi 



dr\ib''/aj \^ -i^Gg '(m^A,7 ' + MiAii^) -mG.vVci"-^ 



+ 



/I -rnGj- 1 (Gi + Mo Ao)/2 \ ( b" /a 



/l + Go'(mn + MiAiGi)/2 J \ib^/a 



+ 



-G^^MoAo - 2AiMo/a mG^^ - + 2m/a)l 

-Go^MiAi 



mGg ^ — vl 



,(18) 



dr \ it 



d h \ _1 {Vg- 4)Go '(m^AQ-^ + M^K^K) -VGGo\m'A.^^ + MiMK) \ f y, 



-m{VG - ^)G^\K + Mo) mVGG^\K + Mo) I \y 



q/Go'MiAi -mGo^\fb"/a\ / 1 + G^ 'MiAiCi/2 -mGoiCo/2 
2 i -mG^^ G^^MoAo J \ib^/a J ^[ -mG^^Ci/2 1 + G^^MoAoCo/2 J \it J ' 



(19) 
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/?=TC. / = 2-^. (22) 



= -AiKy^ - 2AiMoh - 2mit, (20) 
where i/ = 2^/^;, and a; = uo/^GM/K-^ and O = n/^jGM/R^ with M and 7? being the mass and the radius of the star, and 
jr_ d\nMr dlnp F (r/i^)^ 

and 

ci 0)^/31/ ^_ p f_o din a 

B^/8n' ^~ ~dhr7 

Here, equation (16) comes from the continuity equation (3), equations (17) and (18) from the radial and the horizontal 
components of the equation of motion (2), and equations (19) and (20) from the induction equation (4). The matrices Go and 
G\ are defined as 

Go = rnh-o^ - MiAiMq, Gi = m^AJ"^ - MqAoMi, (23) 
and the non-zero elements of the matrices Mq, Mi, Co, and C\ are given by 

W^-.^ = z-Ti-^'T+i' (^o).-..+i = (Ml).-., = |^j;7+i, (Mi),+i,, = |±|j/;+2, (24) 

{C!o)j,j = -{Ij + 2)7;"^+!, iCo)j+i,j = {Ij + 1) {Ci)jj = {Ci)j,j+i = -{Ij + 3) Ji™+2, (25) 

for even modes, and 



(Mo),j = ^Ji7, iMo)j+i,j = j-f^Jr;+i, {Mi)jj = ^Ji;, (Mi),-.,+i = ^J(7+i, (26) 
{Co)j,j = {Ij - i)Ji", (Go)j,j+i = -{Ij + 2) J/'Yi, (Gi)jj = -{Ij + 1) Ji7, (Gi)j+i,j = IjJl'+i, (27) 



for odd modes, and 1 denotes the unit matrix. The definition of the matrices K, Lo, Li, Aq, Ai is the same as that given in 
Lee & Saio (1990). 

The oscillation equations derived above are integrated in the thin surface magnetic layer, and non-magnetic oscillation 
equations (e.g., Lee & Saio 1986) are integrated in the non-magnetic core. The two integrations are matched at the interface 
between the magnetic and non-magnetic regions to obtain a complete solution. 

For numerical computation, we solve a finite set of the linear ordinary differential equations, which is derived by truncating 
the infinite expansions of perturbed quantities so that the first kmax expansion coefficients with Ij and I'j from j = 1 to j = kmax 
are retained for each quantity. In most cases, we employ kmax = 6 in this paper. Using a relaxation method (e.g., Unno et 
al 1989), we solve as an eigenvalue problem for uj the finite set of lineax differential equations with appropriate boundary 
conditions at the center and the surface of the star and jump conditions at the interface (see the next subsection). 



2.2 Bounday conditions and jump conditions 

The surface boundary conditions we use are 

6p/p = 0, ib^ = 0, 6^-1- L+b" = 0, (28) 

where Sp is the Lagrangian perturbation of the pressure, and {L^)ij = Sij{lj + 1). The latter two conditions come from the 
assumption that the perturbed magnetic field must be continuous to the outside field that is given by a scalar potential (i.e., 
V X B' = 0; Cambell & Papaloizou 1986, Cunha & Gough 2000, Saio & Gautschy 2004). The inner boundary condition at 
the stellar center is the regularity condition of the functions and j/j (c-g-, Lee & Saio 1986). 

The jump conditions we apply at the interface are the continuity of the functions j/j and j/2> and the condition that the 
oscillation energy imparted by the perturbed magnetic field is lost from the wave system. To consider the latter condition, it 
is instructive to differentiate equation (18) with respect to Inr to obtain 

d_ d_ C b"/a\ _ a C -G^^MqAo mG];^ - u1\ f -Gq^ Mi Ai rnGg^ \ / b^/a ' 
^'di'-^d^ [ib'^/a ) ^^2 [mG^^^-i^l -Gq^MiAi ) I rnGj"^ -G^HdoAo ) [ib'^/a 



+ ■■■ (29) 



where we have used equations (16) ~ (20) to eliminate h and it. Since the quantity |q| can be very large in the interior, if 
we assume a WKBJ type solution in the region near the interface such that {b^ /a, ib^ /a) oc exp {ikr Inr -|- im4> + iojt), we 
obtain the wavenumber k,. oc ^/a, which suggests the existence of very short magnetic perturbations. These short magnetic 
perturbation are attributable to the Alfven waves. In the deep interior where the gas pressure dominates the magnetic pressure 
so that /3 >> 1, it may be legitimate to write 



ib la j \ib /al„ \ib laj ' \ it I \ it I „ \ it I ' \ y. 
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where the functions with a subscript denote non-Alfvenic perturbations and those with a subscript m denote the Alfvenic 
perturbations, and the former are assumed to have wavelengths much longer than the latter. For the non-Alfvenic perturbations 
in the deep interior, we assume 

mG-' - ul -Gq^MiAi ) \it) ^^^^ \ -,^G-^(rn^Aa^ + M^A-J<) -mG^^/ciQ^ ) \y2 J o'' 
which is equivalent to equations (A4) and (A5) in Loo & Saio (1990), and 

Go'MiAi -mGo^ \ f \ _ d f h\ / 1 + Gq 'MiAiCi/2 -mGo'Co/2 \fh 



(31) 



2\-mGi^ Gi'^MoAoJ yib"^ J ^ dr \it ) ^ \ -mGl^Cxjl 1 + G^'' MoKqCo/2 ) \it 

1 /^(VG-4)Go-i(m%7i+MiAi/s:) -VcG^^Tr? k-^^ + MiAiK)\ f y, 
-m{VG - 4:)G-\K + Mo) mVaG^^K + Mo) J [y 



'2 



(32) 



which comes from equation (4) assuming /3 >> 1 (see §4). Using equations (30) to (32) in equations (18) and (19), we obtain 
for the Alfvenic perturbations 

d_ / b"/a \ _ ( -mG^\Ci + MoAo)/2 \ / b" /a \ f -G-'MqAo mG^^ - z/l \ / h 

\ib'^/a J ^~ \ fl + Go\m'^l + MiAiCi)/2 ) \ib'^/a J \mGo^ - vl -Gq^MiAi ) \it 

f fl -rnG-\C^+ MoAo)/2 \ ( b" /a\ d ( b" /a\ ( b^ /a\ 



r 



fl + Go\m'l + MtAtCi)l2) \ib'^la) ^ dr \ib'^ / a ) ^'^ \ 
^d_(h\ _afG-^MiAi -mGo^\fb"/a\ ^ f 1 + G''' MiAiCi/2 -mGo^Co/2 \fh 

Defining 



dr \it J 2 y -mG;-^ G^^MqAo ) \ib'^ la ) -mGj;'Gi/2 1 + G;"^MoAoGo/2 ) Ut ' ' ^^^^ 




yo = ( ' (35) 

we may formally rewrite equations (33) and (34) as 

r-!^Ym=AYm + ByQ, (36) 
or 

where equations (31) and (32) have been used to eliminate {b^ ,ib^) and (ho,ito), and matrices A and B arc determined by 
equations (33) and (34). Here, we regard equation (36) as a linear differential equation of Ym with an inhomogeneous term 
By^y The homogeneous part of the linear differential equation is solved by assuming Ym oc e'*'''"'', where the wavenumber 
kr is obtained as eigenvalues to the equation ikrYm = AYm- 

To make the discussion more concrete, let us now consider a finite set of the linear differential equations, which is obtained 
by truncating the infinite expansions of perturbed quantities, keeping the first kmax expansion coefficients for each perturbed 
quantity. In this case, the dimension of the vector Ym is 4kmax, and the number of eigenvalues kr to ikrYm = AYm is also 
4:kmax- To construct jump conditions at the interface, we pick up N = 2kmax eigenvalues kl and the eigenfunctions that 
satisfy the condition 

Re(fe^) < and |Re(A;^)| > |Im(fcj)|, (37) 
or 

Im(A;^) < and |Im(A;^)| > |Re(A;^)|, (38) 
and write a general solution to equation (36) as 

Ym = Y^ cjYi,e"^ - A-'By, = {yI, ■■■,yZ)\ [ - A-'By, (39) 

where —A^'^ By^ is a special solution to equation (36), and c-"s are arbitrary constants. The condition (37) means that the 
Alfvenic perturbations are chosen to be running waves into the core (e.g., Saio & Gautschy 2004). It is important to note that, 
different from the case of acoustic modes, the number of the eigenvalues k-j. that satisfy Ro(fcr) < and |Re(fcr)| > |Im(A:r)| is 
usually less than A'^ = 2km.ax for rotational modes. This is the reason why we add the eigenvalues ki that satisfy Im(fcr) < 
and |Im(A;r)| > |Re(fer)| to make the total number of the conditions equal to N. Splitting the solution (39) into the upper and 
the lower parts so that (Ym, • ■ • , ^m)^ and (Ym, ■ ■ ■ , Ym) become square martices of dimension 2kmax x 2kmax, we have 
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Table 1. Complex 
at Q = 0.1. 


eigenvalue 


K = 


ui/Cl and the frequency deviation (kh — kq)/ko of the /' = \m\ r- 


■modes 


Bs{G) KR 


Kl 


m = 


1 m = 2 

{kh - Ko)/ko Xi Ku Kz {kh - Kq) / Kg 




IQi^ 0.99642 
IQi* 0.99642 


2.7 X 10" 
5.0 X 10- 


13 

-9 


1.5 X IQ-" 0.998 0.66517 1.1 x IQ-i" 2.4 x 10"^ 
-1.8 x 10-s 0.98 0.66518 9.5 x lO"'^ 1.8 x IQ-^ 


0.997 
0.98 



{Ym + A-'By,)^ = {Yl,---,YZ)^\ , (40) 

and 

{Ym + A-^Byo)^ = {Yl,---,YZ)A • (41) 

Eliminating the arbitrary constants c' between equations (40) and (41), we finally obtain (e.g., Saio & Gautschy 2004) 

{Yl,---,Y'Z)^' {Ym+A-'By„)^ = {yI, ■ ■ ■ ,Y^)l' {Y,^ + A'' By^) ^ , (42) 

where Ym is to be replaced by Y — Yq with i^o = {bo ,ibo ,ho,ito). The number of jump conditions supplied by equation 
(42) is N, and since the continuity of the functions y^ and j/j interface gives N jump conditions, we have in total 2N 

conditions at the interface. 



3 NUMERICAL RESULTS 

i?-modcs of non- magnetic polytropos were investigated, for example, by Yoshida & Loo (2000a, b). The Coriolis force is the 
restoring force, and the oscillation frequency u) is proportional to the rotation frequency f2. The r-modes are retrograde 
modes, in the corotating frame of the star, having the asymptotic frequency uj = 2mfl/[l'{l' + 1)] for f2 — » for given m and 
and the cigenfunctions are dominated by the toroidal component iT;/ of the displacement vector ^. For isentropic stars, the 
odd r-mode with I' = \m\ is the only r-mode for a given m, and the eigenfunctions are dominated by the nodeless toroidal 
component iT\m\- 

For neutron star models, we employ an A'^ = 1 polytrope of mass M = 1AA4q and radius R = lO^cm, for which we 
assume Fi = 1 + 1/N so that the Schwartzshild discriminant A vanishes in the interior. As discussed in the previous section, 
the radial wavenumber of Alfvenic perturbations near the interface may be approximately given by 

\kr\ ~ ^/a/2 = Uw/2 oc \u>\/Bs, (43) 

where Bs = Bo{R), and tA = v/va with va = \/ B'^/Aivp being the Alfven velocity denotes the traveling time for the Alfven 
waves. This equation indicates that Alfvenic perturbations have very short wavelengths ~ 1/kr if the traveling time of the 
Alfven waves is much longer than the typical oscillation period of the modes. Figure 1 gives the quantity as a function of 
X = r/R for Bs = lO^^G (soUd line), Bs = lO^^G (dashed Une), and Bs = lO^^G (dash-dotted Une), assuming a) = 1. This 
figure shows that, as the depth z = R — r from the surface increases, the radial wavenumber kr of Alfvenic perturbations 
increases quite rapidly and becomes as large as 10'' ~ 10^ beneath the surface, indicating that the travelling time of the Alfven 
waves in a typicl neutron star model with Bs <10^*G is much longer than the oscillation period of rotational modes unless 
~ 0. Since the Alfvenic perturbations associated with the r- modes have very short wavelengths for Bs <10^''G, we have to 
allocate a large number of mesh points to describe them correctly even in the thin surface region. 

In Table 1, we tabulate the complex eigenvalue k = uj/fl = kr + iki and the frequency deviation {kr — ko)/ko for the 
r-modes of m = 1 and m = 2 at Q — 0.1 for Bs ~ lO^'^G and Bs ~ lO^^^G, where kq denotes the eigenvalue of the non- 
magnetic r modes, and Xi = ri/R and r, is the radial distance of the interface from the stellar center. Here, we use kmax = 6 
for Table 1. Note that the location Xi of the interface for the case of Bs = lO^^G is deeper than that for Bs = lO^^G. We 
find that the eigenvalue k_r of the r-modes is hardly affected by the magnetic field of strength Bs <10"'^*G. This is because 
the Alfvenic perturbations become decoupled from the r-modes immediately beneath the surface, giving only a skin effect 
on the modes. If we consider a magnetic field as strong as Bs >10^^, the wavelengths of the Alfvenic perturbations will be 
comparable to the radius of the star for rotational modes of d) ~ f^^l. and no decoupling of the Alfvenic perturbations from 
the rotational modes will be expected in a thin surface layer. In this case, the method of calculation applied in this paper 
may not be justified, but the linearized MHD equations need to be integrated in the entire interior of the star to determine 
the eigensolutions. 
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Table 1 also shows that the r-modes are stable having positive ki, which is, however, quite small for the strength of the 
field considered in this paper. Although the magnetic effects on the eigenfrequency are weak, it is important to note that the 
effects arc larger for the larger Bs, and that the r-mode of m = 1 is less influenced by the dipole field than that of m = 2, as 
indicated by Table 1. 

For the case of Bs = 10^^ G, the real parts of the eigenfunctions of the r-mode of m = 1 are plotted versus r/Rvn Figures 

2 to 4, and those of the r-mode of m = 2 in Figures 5 to 7, whore we assume Cl = Q.l and kmax ~ 6, and we employ the 
amplitude normalization given by Si^ = 1 at the surface. The presence of the magnetic dipole field of strength Bs = lO^^G has 
practically no influence on the modal properties of the modes. The magnetic r-mode has the dominant toroidal component 
iT;/^l„l of the displacement vector over the horizontal _ff;^|^|^i and radial Si^irn\+i components, and the amplitudes of 
the radial and horizontal components are almost the same order of magnitude at the surface. The Alfvenic perturbations 
(6^, ib'^)m, have very short wavelengths and small amplitudes compared with those of {b^ ,ib'^)o- The magnetic perturbations 
(h.it),n themselves have very short wavelengths and very small amplitudes compared to those of {h,ii)a. Wo also find that 
(2/112/2) is hardly affected by the magnetic perturbations near the surface. It may be useful to note an approximate relation 
(h,it)^ ~ (6^,i6^)^ /\/a, which is derived from equation (34) assuming rd/dr ~ \/a for the Alfvenic perturbations. Note 
that in the region near the interface, (6,?^i)o ^constant and (16^2)0 ~ for the r-mode of m = 1, and (6^^2)0 ~ and 
(16^3)0 ^constant for the r-mode of m = 2, where (6^)0 and {ib[)o denote the non- Alfvenic parts of 6^ and ibf , respectively. 
See §4 for a detailed discussion. 

For the case of Bs ~ 10^'* G, the real parts of the eigenfunctions of the r-mode of m = 2 arc plotted versus r/_R in 
Figures 8 to 10, where we assume = 0.1 and kmax = 6, and we employ the amplitude normalization given by Si^ = 1 at the 
surface. The wavelengths and the amplitudes of the Alfvenic perturbations [b^ ,ib'^)m are much longer and much larger than 
those found for the case of Bs = 10^'^G. It is interesting to note that the effect of the magnetic field on the eigenfunctions 
j/j and it becomes apparent near the surface, as shown by the inlets in Figures 8 and 9. In Figure 11, to see the dependence 
of the Alfvenic perturbations on the value of kmax, we plot the real parts of {b^ ,ib^) for the m = 2 r-mode at f2 = 0.1 for 
Bs ~ lO^'^G, assuming kmax = 10. For this case, we have the eigenvalue (k_r — ko)/ko — 1.8 x 10~''' and kj = 8.3 x 10~^, 
which is in good agreement with the value obtained for kmax = 6. Figure 11 shows that the wavelengths and the amplitudes of 
the Alfvenic perturbations {b^ ,ib^)m for kmax = 10 becomes shorter and smaller than those for the case of kmax = 6, which 
may reflect the property of the wavenumber kr determined by ikrYm = AYm that max |Re(A;r)| increases as kmax increases. 



4 DISCUSSIONS AND CONCLUSIONS 

Let us rewrite equation (2) as 



r-Vp' p' „. ^ , 2 (rV xB')xB 

r pgr p r pV Bq 



(44) 



This equation may indicate that in the deep interior where the gas pressure dominates the magnetic pressure such that 
/3 >> 1, the magnetic perturbations B' have practically no influence on non- magnetic oscillations governed by equation (44), 
and they are determined solely by the induction equation (4) . Using equation (4) , we obtain 



Vh ■ B'h = 2Bo (r-^ - l) 
(V, X Bi), = [(,± - 1) (- ™. (V» X 1^)^ - =£) . (. 



cos6>Vif • — -l-sine— + - sin 61— 2 cos 6» — 



dr 



2 



where 



B'h = Bgee + B'^e^, = ^eee -h ^<#.e^, VH=ee-^+e^ 



d 



sin 6 d(j) ' 



Assuming 



for the r-mode of a given m, we obtain from equations (45) and (46) 

H _ 2m ( d \ iTy 
+ \dr ) r ' 



ih. 



21' 
'V + 1' 



d I' 



(45) 
(46) 

(47) 

(48) 

(49) 
(50) 



where J/" = y^P^^rn^)/(4l^ — 1) and I = I' + 1 ior odd modes. Since iT;/ oc r for the r-mode of /' = m = 1, we have 
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{W )o = = constant, {ihi )o — 0, (51) 

and, since iTj' oc for tlie r-mode of Z' = m = 2, we fiave 

(fe^)o=0, {ibJ), = ^l^'IlL^r. (52) 

Equations (51) and (52) for the non-Alfvenic parts of the magnetic perturbations are consistent with the eigenfunctions shown 
in Figures 4 and 7. Note that {bff)o and {ibf)o are both non-zero for the r-modes of /' = m > 2. 

Using a singular perturbation theory, we calculate low-m r-modes of a magnetized neutron star, assuming a magnetic 
dipole field whose axis is aligned with the axis of rotation. For an TV = 1 polytrope of mass — 1.4M0 and radius R — lO^cm, 
we find that the magnetic field of strength Bs <10^*G has no significant influence on the modal properties of the r-modes 
with the frequency lj ~ H unless f2 ~ 0. We show that the magnetic perturbations B' can be divided into Alfvenic and 
non-Alfvenic parts, where the non-Alfvenic part in the deep interior is governed by the induction equation (4). The Alfvenic 
part, on the other hand, has very short wavelengths in the interior. The radial wavenumber kr of the Alfvenic perturbations 
is approximately given by kr ~ tAi^ where tA denotes the traveling time of the Alfven waves, and we have Iai-o >> 1 beneath 
the stellar surface for the r-modes we consider in this paper. The Alfvenic perturbations become very short waves, being 
decoupled from the r-modes in the interior. This decoupling of the Alfvenic perturbations from the r-modes may cause a 
damping effect on the modes. 

If we consider the growthtime t of the r-mode of m = 2 driven unstable by current multipole gravitational radiation, we 

have 

r ~ {3/4fn=2^~'^ ~ 1.4^"® sec, (53) 



where we have used fi=2 = 3.31sec, and Q = Q/ao with ctq — ^ GMjB? (see, e.g, Yoshida & Lee 2000). Defining the growth 
timescale as r = —l/u>i, we have 

K/ ~ -5 X 10~^^^ (54) 

where we have used cro = 1.36 x 10*s"^ for M = 1.4M0 and R = 10*^ cm. Thus, we have ~ -5 x 10"^" for Cl = 0.1 and 
Ki ~ —1.2 X 10~^ for Cl — 3/4. Comparing these values with the numbers ni given in Table 1, we may find that the r-mode 
instability can be damped by dissipative very short Alfvenic perturbations in the surface layer when the neutron stars with 
Bs » 10^'^G are slowly rotating. For neutron stars in LMXB's, however, the strength of magnetic field is believed to be 
much weaker than lO^^G, and the effects of the magnetic field on the r-mode instability will be unimportant. 

In this paper, we have assumed that the axis of the magnetic field is aligned with the axis of rotation. If we lift this 
assumption, the problem will be much more complicated since the axisymmetry cannot be assumed and hence a perturbed 
quantity must be represented by a sum of terms proportional to Yj™ over the indices I and m. The possible existence of a 
crustal layer in the envelope of a rotating magnetized neutron star also makes difficult the problem of the oscillations (e.g., 
McDermott, Van Horn, & Hansen 1988; Lee & Strohmayer 1996). Magnetic effects on the torsional oscillations of a neutron 
star with a solid crust have been discussed by Carroll et al (1986) and Messios, Papadopoulos, & Stergioulas (2001), but no 
effects of rotation have been included in these studies. If there exists a fluid ocean on the solid crust of a rotating neutron 
star, the magnetic effects on waves propagating in the ocean may be important observationally if the waves are responsible 
for burst oscillations found in neutron stars in LMXB systems (e.g., Lee 2004). 

Acknowledgments: The author is very grateful to Prof. H. Saio for his illuminating discussion on the problem. 
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r/R 

Figure 1. Plot of a^^^ as a function of r/R for an A'^ = 1 polytrope of mass M = IAMq and radius R = lO^cm, where w 
is assumed. Solid, dashed, and dash-dotted lines are for the cases of Bs = IQi^G, 10"G, and lO^^G, respectively. 




r/R 

Figure 2. Eigcnfunctions (''yi);=|m|+i (solid curve) and {i'y2)i=\m\+i (dashed curve) as a function of r/R for the r-mode of m = 1 
CI = 0.1, where Bs = lO^^G and kmax = 6. 
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r/R 

Figure 3. Eigenfunctions iTj/_|^| (solid line) and ^^;=|to|+i (dashed line) as a function of r/R for the r-mode of m = 1 at H = 0.1, 
where Bg = lO^^G and kmax = 6. 
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Figure 4. Eigenfunctions (lower curve) and (upper curve) as a function of r/R for the r-mode of m = 1 at Q = 0.1, 

where Bs = 10^^ G and kmax = 6. 




Figure 6. Same as Figure 3 but for m = 2. 
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Figure 7. Same as Figure 4 but for m = 2. 
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Figure 8. Same as Figure 5 but for Bs = lO^^G. 
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Figure 11. Same as Figure 10 but for kmax = 10. 



